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ABSTRACT 

We revisit in more detail a model for element abundance fractionation in the 
solar chromosphere, that gives rise to the "FIP Effect" in the solar corona and 
wind. Elements with first ionization potential below about 10 eV, i.e. those that 
are predominantly ionized in the chromosphere, are enriched in the corona by a 
factor 3-4. We model the propagation of Alfven waves through the chromosphere 
using a non-WKB treatment, and evaluate the ponderomotive force associated 
with these waves. Under solar conditions, this is generally pointed upwards in 
the chromosphere, and enhances the abundance of chromospheric ions in the 
corona. Our new approach captures the essentials of the solar coronal abundance 
anomalies, including the depletion of He relative to H, and also the putative 
depletion of Ne, recently discussed in the literature. We also argue that the FIP 
effect provides the strongest evidence to date for energy fluxes of Alfven waves 
sufficient to heat the corona. However it appears that these waves must also 
be generated in the corona, in order to preserve the rather regular fractionation 
pattern without strong variations from loop to loop observed in the solar corona 
and slow speed solar wind. 

Subject headings: Sun: abundances - Sunxhromosphere - turbulence - waves 



1. Introduction 

Since about 1985 it has been recognized that the composition of the solar corona and 
the photosphere are not the same. In the corona, the ratios of abundances of elements with 
first ionization potentials (FIP) less than about 10 eV relative to abundances of elements 
with a FIP greater than about 10 eV are about a factor of 3-4 times higher than in the 
photosphere. Those elements with a FIP greater than 10 eV appear to have a photospheric 
composition in general (with respect to hydrogen) in the corona, and the low FIP elements 



are enhanced in abundance there. This fractionation has recently been explained ( ILaming 
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2004al ) as being due to the ponderomotive force in the chromosphere from Alfven waves. This 
is usually directed upwards, and acts on chromospheric ions, but not neutrals. Elements that 
are predominantly ionized in the chromosphere (low FIP elements like Al, Mg, Si, Ca, and 
Fe) can be enhanced in abundance as they flow up to the corona, whereas high FIP elements 
such as C, N, O, Ne, and Ar that are largely neutral appear essentially unaffected. The 
abundance of sulfur (FIP = 10.36 eV) is between photospheric and coronal. 

It was recognized very early that the most plausible site for FIP fractionation was the 
chromosphere, where low FIP elements are generally ionized and where high FIP elements 
are at least partially neutral. Henoux (1995, 1998) reviewed the early models. Those which 
rely on ion-neutral separation by diffusion in magnetic fields, in temperature gradients, or in 
upward plasma flows suffer from problems of the speed of the process (diffusion is inherently 
slow) or the choice of boundary conditions. Realistic FIP effect models must include some 
form of external force that acts upon the plasma ions but not upon the neutrals. The first 
effort along these lines (Antiochos 1994) considered the cross-B thermoelectric electric field 
associated with the downward heat flux carried by electrons which gives rise to chromospheric 
evaporation. This draws ions into the flux tube, enhancing their abundance over neutrals. 
The absence of a FIP effect in coronal holes arises naturally, but in coronal regions where 
FIP fractionation occurs, a mass dependence is predicted, which is not observed. 

Henoux & Somov (1992) proposed that cross-B pressure gradients arising in current 
carrying loops could enhance the ion abundance by a "pinch" force. Azimuthal motions 
of the partially ionized photosphere at flux tube boundaries generate a system of currents 
flowing in opposite directions, such that the azimuthal component of the field vanishes at 
infinity. Details of the fractionation (mass dependence, degree, etc) remain to be worked out, 
but it is thought to begin just above the temperature minimum region at about 4000 K, and 
continue until temperatures where all elements are ionized. More recent suggestions have 
been that chromospheric ions, but not neutrals, are heated by either reconnection events 
(Arge & Mullan 1998), or by waves that can penetrate down to loop footpoints from the 
corona (Schwadron, Fisk & Zurbuchen 1999). Strengths and weaknesses of these two models 
are discussed in some detail by Laming (2004). For the time being we comment that although 
both models seem capable of producing mass independent fractionation of about the right 
degree, they, in common with the others mentioned above, only predict positive FIP effects. 
In these models there is no possibility of an "inverse FIP" effect such as seen in the coronae 
of active stars (see e.g. Feldman & Laming 2000, Laming 2004 and references therein). 

This consideration led Laming (2004) to consider the action of ponderomotive forces due 
to Alfven waves propagating up through the chromosphere and either transmitting into or 
being reflected from coronal loops. This force can in principle be either upward or downward 
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and is given approximately by (see derivation in Appendix A) 

F _ q 2 dSEj 
4m (Q 2 — uj 2 ) dz 



X 



where Q and u are the particle cyclotron frequency and the wave frequency respectively, q 
and m the charge and mass, and E± is the wave peak transverse electric field. The depen- 
dence on the Alfven speed, Va, means that the ponderomotive force is usually strongest at 
the top of the chromosphere. The ponderomotive acceleration, F/m, is independent of the 
ion mass, leading to the essentially mass independent fractionation that is observed. Using 
an analytic model solar loop from Hollweg (1984), upward ponderomotive forces on ions are 
much more common in solar conditions, and for typical density gradients in the chromosphere 
can be larger than the gravitational force downward. The magnitude of the FIP fractiona- 
tion is dictated by the resonant properties of the coronal loop and the corresponding wave 
energy density. Large loops with resonant frequencies similar to the chromospheric period of 
200-300 seconds admit strong Alfven wave fluxes with correspondingly large FIP fractiona- 
tion. Open field lines, or coronal holes, with formally infinite wave periods or unresolved fine 
structures (i.e. the chromosphere and lower transition region away from coronal loop foot- 
points, Feldman 1983, 1987) with much shorter wave periods do not have much Alfven wave 
transmission into them and so would have very low FIP fractionation, as observed. Using 
model chromospheres from Vernazza, Avrett, & Loeser (1981), most of the FIP fractionation 
is found to occur at the top of the chromosphere at altitudes with strong density gradients 
near the plateau regions, where low FIP elements are essentially completely ionized and high 
FIP elements are typically at least 50% ionized. The Laming (2004) model thus comes about 
as a natural extension of existing work on Alfven wave propagation in the solar atmosphere 
with essentially no extra physics required. 



In this paper we revisit the lLamingi (j2004al ) model using a numerical treatment of Alfven 



wave propagation in a coronal loop rooted in th e chromo s phere at each footpoint. Section 2 



describes this model and the improvements over iLamingi (l2004al ) with a series of illustrative 
examples. Section 3 gives a more complete tabulation of the FIP fractionation expected in 
a variety of elements, and section 4 gives some discussion of the implications of the new 
models, both for the solar coronal abundances of helium and neon, and for MHD wave origin 
and propagation in the solar atmosphere, before section 5 concludes. 
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2. Ponderomotive Driving of the FIP Effect 

2.1. Introduction and Formalism 

Laming (2004) used a WKB approximation to treat the case of strong transmission of 
Alfven waves into a coronal loop, and hence evaluate the FIP fractionation. Here we have 
extend this work to a no n- WKB treatment. The procedure follows that described in detail 
by Cranmer & van Ballegooijen (2005), but applied to closed rather than open magnetic 
field structures. The transport equations are (see Appendix B for a derivation) 

^ + ( „ ± ^ (m±w (_|_ + _|_), (2) 

where I± = Sv ± SB/^JAitp are the Elsasser variables for inward and outward propagating 
Alfven waves respectively. The Alfven speed is Va, the upward flow speed in u and the 
density is p. The signed scale heights are H D = pj (dp/dz) and Ha = Va/ (OVa/Oz). In the 
solar chromosphere and in closed loops we may take u « Va- 

We use the same chromospheric model as before in Laming (2004), but this time als o 



embedding it in a 2D force- free magnetic field computed from formulae given in lAthayl (119811 ). 
and shown in Figured! We take a scale here of 1 unit to 1,000 km, and place the bottom of 
the plot at 500km altitude in the chromosphere. The region where the plasma beta (ratio of 
gas pressure to magnetic pressure) equals unity is taken at 650 km altitude, or at y = 0.15 
in the figure. This region is where upcoming acoustic waves from the convection zone can 
convert to Alfven waves by mode or parametric conversion, or where downgoing Alfven 
waves can convert back to acoustic waves. We adopt an altitude of 650 km as the boundary 
of our simulations where Alfven waves are launched upwards. We consider a loop model 



similar to that of iHollwegi (1l984j ). the coronal portion of which is illustrated in Figure [2j 
Waves from beneath impinge on the right hand side chromosphere-corona boundary and are 
either reflected back down, or transmitted into the loop. Waves in the loop section bounce 
back and forth, with a small probability of leaking back into the chromosphere at either 
end. The magnetic field in the coronal loop section is uniform, and is simply the extension 
of chromospheric force-free field into the corona. In the current work, the non-uniformity 
of the magnetic field in the low chromosphere has no effect on the FIP effect, since the 
fractionation in this work appears towards the top of the chromosphere where the magnetic 
field is almost parallel. The loop density is similarly extrapolated, but with a density scale 
height taken here to be equal to the loop length. All models presented here assume waves 
propagating up the y-axis at x — 0. 

Equations (2) are integrated from a starting point in the left hand side chromosphere 
(hereafter chromosphere "A") where Alfv'en waves leak down into the chromosphere, back 
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through the corona to the right hand side (chromosphere "B") where waves are fed up from 
below. In this way the reflection and transmission of Alfven waves at the loop footpoints 
and elsewhere is naturally reconstructed. The velocity and magnetic field perturbations are 
calculated from 
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The wave energy density and positive and negative going energy fluxes are 
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and the wave peak electric field appearing in equation (1) is 
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Given the ponderomotive a cceleration from equation (1), the FIP fractionations are 
calculated in a similar manner to iLamingl (2004a), but wi th one important modification. In 
the momentum equations (10) and (11) of lLamingi (|2004al ). we add the motion in the wave to 



the ion and neutral partial pressures, so that P s ,i, n = {k&T/m s + vf urb + v%, ave ) p s ^ n /2 where 
the first two terms in parentheses represent the ion thermal velocity and the microturbulent 
velocity in the model chromosphere, and the third term represents the motion of the ion in the 
Alfven wave. In true collisionless plasma, neutrals would not respond to the wave. However 
the solar chromos phere is sufficiently collisional that neutrals move with the ions in the 
wave motion (e.g. IVranjes et al.ll2008l ) for wave frequencies well below the charge exchange 
rate that couples neutrals and ions, and so neutrals require the same form for their partial 
pressure as the ions. Followin g the derivation through, writing u s>n dP S) i/dz + v si dP Stn jdz 
from equations (10) and (11) in lLamingi (l2004al ) (with the ponderomotive term in dp s ^/dz 
0, i.e. b = 0) we find 



d_ 

dz 
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+ V turb + v wi 



+ p s [g + u e ff (u s - u)] + p s a£ s v e ff/v si = 



(6) 



where u eff = v s ^v s ^j %v s , n + (1 - £ s 



fx 



with £ s being the ionization fraction of element 



s, and v s ^ and v s ^ n the collision rates of ions and neutrals respectively of element s with the 
ambient gas, and a the ponderomotive acceleration. This leads to 
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where v 2 s = kT/m s + v^ urb + v^ ave . We argue that in the absence of the ponderomotive 
acceleration a, the effect of turbulence should be to fully mix the element composition. Thus 
we choose u s in equation (7) such that 



V s {Zif 



-9 - "eff {Us ~ U) 



dz 



to yield the fractionation by the ponderomotive force as 

Ps (gu) _ ^ f f Zu Zsais ef f/vs,i 



dz 



(9) 



As mentioned previously (jLamingl l2004al ) . the solar chromosphere is undoubtedly more 
active and dynamic than represented by equations (6-9). However this choice allows us to 
model a chromosphere which in the absence of Alfven waves is completely mixed, presumably 
by hydrodynamic turbulence, and upon which the ponderomotive force acts to selectively 
accelerate chromospheric ions. Chromospheric simulations excluding turbulence find huge 
(and unobserved) variations in various corona l element abundances due to ambipolar and 
thermal diffusion (IKillie fc Lie-Svendsonl 120071 ). It would appear that such abundance vari- 
ations are unavoidable in the situation that the ch romosphere remains undist urbed for a 
sufficient length of time (days to weeks in the case of IKillie fc Lie-Svendsonl 120071 ). We argue 
that hydrodynamic turbulence acts on timescales much shorter than this (but still longer 
than that required to establish the FIP effect) to leave a fully mixed chromosphere in the 
absence of ponderomotive forces. We model the ponderomotive acceleration as acting only 
on the chromospheric ions, since the ponderomotive acceleration divided by the flow velocity 
a/u ~ 10 — 1000 s^ 1 is much greater than the charge exchange rate (of order 1 s _1 ). This 
rate should also be sufficiently greater than the turbulent mixing rate discussed above. 



We also use a more recent chromospheric mod el (model C7 in 



introduced as an update to t he older VALC model (jVernazza. Avrett. &: Loeserlll98ll ) used 



m 



Avrett fc Loesei 



2008|) 



previously in iLaminel J2004ah. We continue to evaluate photoionization rates using incident 
sp ectra ba s ed on Vernazza fc Reeved (119781 ) with the extensions and modifications outlined 



Laming! (j2004al ). In most cases, the "very active region" spectrum is used, being the most 



consistent with the underlying chromospheric model. 



2.2. A Loop On Resonance 



In the following, Figures [31 HI and [5] illustrate the solutions for a loop 100,000 km long, 
with coronal magnetic field B = 9.9 G. This yields a wavelength for 3 minute period waves 
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approximately the same as twice the loop length, and therefore such waves can be transmitted 
into the coronal loop section from the chromosphere. We concentrate on 3 minute waves 
since u nlike 5 minutes waves, these require no special conditions to propagate up into the 
corona (Ide Pontieu et al.l 120051 ). Figure [3] shows the coronal section of the loop. From top 
to bottom the three panels give the amplitudes of 5v and dB/y/Anp in units of km s _1 . Real 
and imaginary parts are given as black and gray lines respectively, with 5B/y/47rp and Sv 
given by solid and dashed line respectively. The wave amplitude has been chosen to be ~ 30 
km s" 1 in the corona, giving a typical spectral line FWHM consistent with observations (e.g. 



Mcintosh et al.ll2008l ). The second panel gives the wave actions (energy fluxes) for the left 



going (solid) and right going (dashed) lines, and their difference divided by the magnetic field 
as a dotted line, in arbitrary units. This last quantity should be a straight horizontal line 
if energy is properly conserved in the calculation. The third panel gives the ponderomotive 
acceleration, in cm s -2 . Throughout the coronal section of the loop, it is significantly lower 
than the gravitational acceleration. Solid lines indicate positive, i.e. right going, and dashed 
lines indicate left going accelerations. The osc illation amplitud e has been chosen to give 
mass motions within observational constraints (IChae et al.lll998l ; iMcIntosh et al.ll2008l ). as 
measured from line profiles. 

Figure H] shows the same three plots for the left hand side chromosphere "A" , where 
waves leak down from the corona, together with a fourth panel showing the degree of frac- 
tionation for the abundance ratios Fe/H, O/H, and He/H. The ponderomotive acceleration 
in the chromosphere is much larger than in the corona, especially towards the top. The 
most significant fractionation occurs here, increasing Fe/H in this case by a factor of 1.4 over 
photospheric values, O/H by a factor of around 1.25, with He/H remaining nearly unchanged. 

Finally Figure [5] shows the same four panels for the chromosphere "B" on the right 
hand side, where the upgoing Alfven waves originate. The ponderomotive acceleration is 
still pointed upwards (though is negative in the coordinate system used here), giving the 
same FIP fractionations as before. In exact resonance, the chromospheric ponderomotive 
force behaves the same as at the opposite footpoint already shown in Figure HI 



2.3. A Loop Off Resonance 

Figures El CD and M show the same variables as before, but for a loop 100,000 km long 
and with magnetic field B = 19.8 G. Now the loop is a quarter wavelength long, and almost 
complete reflection of the incident Alfven waves on the right hand side takes place. The 
simulation has been normalized so that the incident Alfven wave flux coming up from the 
chromosphere is about the same as for the on resonance case. Figure M shows that the 
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coronal loop oscillation is now much weaker than before, by about a factor 20. In the left 
hand chromosphere "A" (Figure [7j) negligible Alfven wave flux leaks through and no FIP 
fractionation occurs. In the right hand chromosphere "B" (Figure ED, the behavior is quite 
different to the previous case. The ponderomotive force is now downwards pointing for 
most heights in the chromosphere, and it is still very small, also giving essentially no FIP 
fractionation. The downward directed ponderom otive for c e migh t be of interest in cases 
where the turbulence is stronger. As reviewed in iLamingi (j2004al ). the coronae of various 
active stars exhibit an inverse FIP effect, where the low FIPs are depleted in the corona 
instead of being enhanced. The reversal of the ponderomotive force under these conditions 
is a plausible mechanism for such abundance anomalies. 



2.4. Loops with Stronger Turbulence 

The first case above was designed to give a coronal nonthermal mass motion within 
observational limits, i.e. a root mean square 8v ~ 30 km s _1 . The upgoing energy flux 
of Alfven waves at the loop footpoint is ~ 10 5 ergs cm -2 s -1 , and is insufficient to power 
the coronal radiation power loss by one to two orders of magnitude. In this subsection we 
consider the same loop as in the first case, but with an Alfven wave upward energy flux of 
about 2 x 10 6 ergs cm -2 s -1 ; sufficient to power radiation from a 100,000 km loop with a 
density of 10 8 — 10 9 cm~ 3 . The predicted nonthermal mass motions in the corona are now 
unphysically high, in excess of 100 km s _1 , unless we are able to argue that only a small region 
of the corona oscillates with this speed. We discuss this further in subsection 4.1. This is less 
of a problem in the transition region where the "classical" transition region that conn ects a 
coron al loop with the chromosphere has only recently been ide ntified in observations (Peter 



20011) . being otherwise masked by "unresolved fine structures" (jFeldmanlll983l . 119871 ). 



Peter 



(120011 ) in fact observed nonthermal line broadening in what he interprets as the "classical" 
transition region approaching the values modeled in this section, and suggests that they arise 
from the passage of an Alfven wave with sufficient energy flux to heat the corona. In this 
strong wave field, the behavior of the FIP fractionation is now subtly different, as shown in 
Figures 191 and [TOl for the left and right hand side chromospheres ("A" and "B") respectively. 
On each side, Fe is somewhat more enhanced than in the previous case, at 3 - 3.3. O has 
a FIP fractionation of about 1.7-1.8, similar to before, but He/H is now at about 0.8 of its 
photospheric value. 

This new behavior can be understood with reference to equation 5, and the denominator 
in the integral, v 2 s = kT/m s + vf urb + f 2 a „ e . In the case that the first two terms in t> 2 
dominate, i.e. weak Alfven turbulence, all elements (high FIPs as well as low FIPs) are 
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fractionated positive to H because H has the largest thermal velocity in the denominator. 
When the Alfvenic velocity dominates, the fractionation changes and is determined solely 
by the numerator in the integral. In this case the element that stays neutral the longest, He, 
as expected, has the lowest abundance in the corona, being depleted with respect to H. This 
occurs because H experiences a stronger ponderomotive enhancement. O/H is unchanged, 
again as expected because O and H have very similar ionization potentials and their ionization 
structures are locked by charge exchange reactions between them. Fe remains fractionated 
with respect to H, by a similar amount as before. The inclusion of the Alfven turbulence in 
v s leads to a natural saturation of the FIP effect, at about the level observed. Thus for a 
wide range of turbulence levels, and FIP effect of around 3 should be expected. 

The decrease in He/H is especially int eresting. It might b e relevant to the He abundance 



in the solar wind, of around 4-5% (e.g. lAellig et al.l l200ll ; iKasper et al.l 120071 ) compared 
with a photospheric abundance of 8%, also seen in coronal holes and quiet solar corona 
( jLaming fe Feldmanl 1200 ll . 120031 ). and is discussed further below. 



3. More Realistic Examples 



We put three Alfven waves with angular frequencies 0.025, 0.022, and 0.016 rad s _1 , with 
relative intensities 1:0.5:0.25 in the left h and chromosphe re designed to match the network 
power spectrum displayed in Figure 1 of iMuglachl (120031 ). The loop is 100,000 km long as 
before, with a magnetic field of 7.1 G, which puts the 0.025 rad s _1 on resonance. In the first 
case, FIP fractionations a re computed for the left h and chromosphere "A", using the very 
active region spectrum of IVernazza fe Reeved (119781 ). and are given in Table 1. This is the 
region where waves leak into the chromosphere from the corona before being reflected back 
up again, and should give FIP fractionation. The corresponding model is shown in Figure 
11. There are three important differences from the tabulation given previously (Table 2 in 
Laming! l2004al ). The first is that with increasing wave energy flux, the FIP fractionations 
now appear to saturate at levels corresponding to a fractionation of low FIP elements over- 
abundant with respect to high FIP elements by a factor of about 3, and does not increase 
without limit. This arises from the inclusion of the term in v\ ave in the ion and neutral par- 
tial pressures discussed above, and means that for a wide range of turbulent energy densities, 
similar fractionated abundances should result. The other new features, already mentioned 
briefly above, are the depletion in the He abundance, and at higher energy fluxes also the 
Ne abundance relative to H. These also stem from the modification to the partial pressures. 



Th ese new calcul a tions are co mpared in Table 1 with ob ser vations fromlZurbuchen et al. 



(120021 ). iBryans et al.l (120081 ) and iGiammanco et al.l (120081 ). IZurbuchen et al 



(12002( 1 give 
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abundances measured in the slow sp eed solar wind during 19 97 /8 relative to O, re lative 
to photospheric abundances given by iGrevesse fe Sauvall (119981 ). iBryans et al.l (120081 ) give 
abundances observed spectroscopically in a region of quiet sol a r coro na, again tabulated 
relative to the photospheric composi tion of IGrevesse fe Sauvall (119981 ). with small modifi- 
cations by iFeldman fe Laming! (120001 ) . With the exceptions of Mg and K, the calculated 
abundances agree well with those observed for a wave energy flux between one and four 
times that shown in Figure 11, both for the elements that are depleted, like He and Ne, and 
for those enriche d. We p r edict s tronger fractionation in Mg than is in fact observed, and 
stronger than in iLamingl (l2004al ). The reason for this h as been tracked to the use of the 
newer chromospheric model from Avrett fe Loeser (|2008 ). where H retains a higher degree 
of ionization lower in the chromosphere than in the previous VAL models. This then in turn 
renders the ionization of Mg probably spuriously high because of charge transfer ionization 
with the ambient protons. The difference in ionization fraction between 0.99 and 0.95 makes 
a considerable impact on the fractionations that result. Other low FTP elements, Si, and Fe, 
do not have charge transfer ionization rates tabulated by lKingdon fe Ferlandl (119961 ). and so 
are unaffected by this change. The cause of the discrepancy for K is less clear. Like Na, K 
is very highly ionized throughout the c hromosphere due to its very low FIP, and should be 
expected to fractionate strongly, though IBryans et al.l (120081 ) do comment that their analysis 
only includes one line of K IX. 

The results of the calculation for the right hand side chromosphere "B" are given in 
Figure 12. The loop model chosen is resonant with the 0.025 angular frequency wave, 
and this is the component transmitted into the corona. However the FIP fractionation is 
significantly reduced by the presence of the other wave frequencies which are reflected from 
the corona. In the right hand chromosphere "B" the weaker components on the left are now 
the strongest. This does not produce much change in the ponderomotive acceleration, but 
increases the term u^ ai)e in the denominator of the integrand in equation 6, thereby reducing 
the fractionation. A wave source in the chromosphere is unlikely to be monochromatic, and so 
this situation of partial transmission and partial reflection with the reduced FIP fractionation 
will be ubiquitous in the solar atmosphere. This does not agree with observations, for which 
chromosphere "A" is a much better match. We therefore argue that if the FIP effect is due 
to the ponderomotive force of Alfven waves in the chromosphere, these must have a source 
in the corona. We return to this thought in subsection 4.1. 

Although this calculation has been done for a closed loop, we expect that this chromo- 
spheric wave pattern will also arise at the footpoint of an open field line in a coronal hole. In 
fact the character of our chromospher i c solu tion matches well with t hat found in the o pen fi eld 
case by lCranmer fe van Ballegooijenl (120051 ). subsequently shown in lCranmer et al.l (120071 ) to 
exhibit FIP fractionation similar to that observed in the fast solar wind. We emphasize this 
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point by showing in Figures 14 and 15 the coronal and chromospheric portions of an open 
field flux tube. In this case we start the integration at an altitude of 5 x 10 5 km with purely 
outgoing waves, and work back to the solar surface. This restricts us to the region where 
the solar wind outflow speed is still much lower than the Alfven speed, in keeping w ith our 

Jl998h . mod- 



assump tion of u « Va above. We take mag netic field from lBanaskiewicz et al. 



ified by lCranmer fc van Ballegooijenl (120051). and choose a density scale height to match the 
observed and modeled density profiles in lLaming ( 2004b ). Figu re 14 shows Sv and SB/ \/4ttp, 
chosen to match the observational and modeling constraints in ICranmer fc van Ballegooijen 
(120051 ). Figure 15 shows the extension of these variables into chromosphere "B", in a sim- 
ilar manner to the previous figures. While there is much more to be said about the wave 
properties in open field lines, the important point to be made here is that the ponderomo- 
tive force naturally produces a very small fractionat ion in this geometry. T his is consistent 
with o bserved abundances in the fast solar wind ( Zurbuchen et all I2OO2I ) and in coronal 
holes (IFeldman et al.l 119981 ). A tabulation of coronal hole fractionations is given in Ta- 
ble 2, in a similar f ormat to that in Table 1, using the coronal hole incident spectrum of 



Vernazza fc Reeved (119781 ). 



4. Discussion 



4.1. Alfven Wave Energy Fluxes 



It appears from the forgoing that the abundance anomalies observed in various regions 
of the solar corona may yield inferences on the energy fluxes of Alfven waves in the chro- 
mosphere. Our initial considerations then imply that wave energy fluxes sufficient to heat 
the solar corona or accelerate the solar wind are necessary to produce the correct fraction- 
ation. Energy flu xes observed in slow m ode and fast mode waves are not sufficient to heat 
the solar corona (lErdelyi fc Fedunl 120071 ). The detection of Alfven waves, the favored mode 
for transporting energy to the solar corona, is much harder, since they are incompressible 
and can only be revealed through Doppler shifts or motions, which become hard to see in 
inhomogeneous conditions where Alfven waves on nei ghboring flux su r faces can propagate 
at different speeds and lose phase coherence. However iTomczyk et al.l (120071 ) claim the de- 
tection of Alfven waves in the s olar c orona, albeit with insufficient energy flux to heat the 
corona. Ivan Doorsselaere et al.l (120081 ) argue that the detected waves are in fact kink mode 



waves, for the reasons suggested above. Ide Pontieu et al.l (120071 ) observe transverse waves in 
the chromosphere, which they argue should be interpreted as Alfven waves in the absence of 
a chromospheric waveguide. However these waves are inferred from the observed oscillations 
of spicules, which clearly have radial structure. The energy flux detected by these authors 
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~ 10 5 ergs cm 2 s 1 is close to being sufficient to heat the solar corona or accelerate the 
solar wind. 

In this paper we argue that the FIP effect is due to the ponderomotive force associated 
with transverse waves in the chromosphere. Longitudinal MHD waves do not generate electric 
field. In order to generate the observed FIP fractionation, the energy fluxes associated with 
these waves need to be of order 10 6 — 10 ergs cm -2 s _1 , much closer to those required for 
coronal heating. It is clear that FIP fractionation is associated with the transmission of waves 
between the chromosphere and the corona, and correspondingly we argue that the required 
transverse waves should be identified as Alfven waves to meet this condition. The fast mode 
totally internally reflects somewhere in the transition region or low corona (Schwartz & Leroy 
1982, Leroy & Schwartz 1982). 

The nonthermal mass motions predicted in the coronal section of this loop are higher 
than observed. In the transition region, this is not necessarily a problem since ample evidence 
exists to show that the "classical" transition regions of coronal loops are rarely observed. 



1983. 


1987; 


Peter 


2001) 



filament or shells as in Alfven resonance models (e.g. iTerradas et al.l 120081 ) while the rest 
of the emitting loop undergoes much slower oscillations. Another possibility might be the 
generation of turbulen ce following coronal reconnection events associated with nanoflares 
(IDahlburg et al.ll2005l ). In each case it is likely that the turbulence would actually be pro- 
duced in the coronal section of the loop, not in the chromosphere, and will also be in 
resonance with the loop. This also appears to be the conclusion to be drawn from section 3. 
Chromosphere "A" where waves leak down from the corona before being reflected back again 
gives stronger and more consistent FIP fractionations for a wide variety of wave spectra than 
chromosphere "B", where waves are incident upwards on the loop from the chromosphere 
below. We therefore argue that the FIP effect is more likely to arise with a coronal source o f 
Alfven waves, rather than a chromospheric source as originally conceived in lLamingi (j2004al ). 



and that this inference will constrain the means by which the corona may be heated. 



4.2. Helium and Neon in the Solar Corona 



The fractionations computed in this paper differ from those in iLamingi (j2004al ) in three 
notable ways. First, as the turbulent energy density increases, the fractionation does not 
increase without limit but saturates at values broadly consistent with those observed. This 
is due to a refinement in our formalism discussed above, where the wave oscillation velocity 
is included in the ion and neutral partial pressures in equation 3. 
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The second is that at high turbulence levels, He becomes significantly depleted relative to 
H. Comparing Tables 1 and 2, we find a stronger depletion in the coronal loop, representative 
of the slow speed solar wind, than we would in a coronal hole, the source of t he fast wind. 



The abundance ratio He/H in the fast solar wind is fairly constant at about 5% (jAellig et al. 



200ll ; IZurbuchen et al.l l2002f ) . or a depletion of 0.59 from the phot ospheric val u e of 8 .5%. 
He/H in the s low sp eed solar is lower, and generally more variable. lAellig et al.l (120011 ) and 
Kasper et al.l (120071 ) find He/H varying with wind speed, with these variations being more 
pronounced at solar minimum, where He/H ~ 1% for speeds below 300 km s^ 1 , approaching 
4.5% for spee d s abo ve 500 km s _1 . At solar maximum, He/H is always in the range 3.5 - 5%. 
Kasper et al.l (120071 ) also find a dependence on heliographic latitude during periods of solar 
minimum, with lower He/H being found closer to the heliographic equator. Table 1 gives 
values of He/H down to about 3.5%. Overall though, our modeled values for the abundance 
ratio He/H are very encouragingly consistent with observations, lending confidence to our 
approach. 

The third, and most contro yersial new feature is the similar depletion predicted for Ne. 
This was originally suggested by lDrake fc Testa! (120051 ) from a survey of the Ne/O abundance 
ratio in a sample of late-type stellar coronae, as a solution to the problem in helioseismology 



presented by the reduction in the sola r photospheric abun dance of O (see lCaffau et al.ll2008 



and references therein). Specifically, (IBasu fc Antial 12004 ) the depth of the solar convection 
zone demands a metallicity higher than that coming from the st andard solar compos ition, 
with the O abundance rey ised downwards by nearly a factor of 1.5 (lAsplund et al.ll2004l ) from 
Grevesse fc Sauvall (119981 ) . Ne, having no photospheric lines on which to base an abundance 
measurement, was suggeste d as the element most likely to resolve t his by having a higher 
postulated abundance (e.g. iBahcall et al.ll2005l ; IBasu fc Antial 120081 ) . iDrake fc Testal (120051 ) 
find coronal Ne/O typically ~ 0.4 in stars which exhibit either no FIP effect or an inverse 
FIP effect, and argued that the general consistency of Ne/O among their sample of 21 stars 
suggests no significant fractionation between Ne and O here between pho tosphere and corona . 



The solar co ronal abundance ratio Ne/O, measured at ~ 0.15 — 0.18 (jSchmelz et al.l 12005 



Young! 120051 ) would imply therefore that Ne is depleted in the solar corona relative to the 
photosphere, similarly to He. Our calculations in Table 1 provide some support to this view, 
especially at higher turbulence levels, where Ne/O is about 0.5 of its photospheric value. 



4.3. Fractionation in the Low Chromosphere 



One main feature of Laming (2004) model and the calculations presented above is that 
the fractionation is predicted to occur relatively high up in the chromosphere, at altitudes 
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greater than 2000 km. However in the literature there are already indications that, at least 
in active regions and flares, that fractionation should set in lower down. 

In an analysis of HRTS II (the Naval Research Laboratory's High Resolution Telescope 
and Spectrograph) data, Athay (1994) observed variations in the C I 1561A /Fe II 1563 A 
line intensity ratio. Compared to plage regions around a sunspot, the sunspot itself has a 
higher ratio C I/Fe II, while surrounding C I dark flocculi have a lower ratio. Similar results 
are found by Doschek, Dere & Lund (1991) and Feldman, Widing & Lund (1990). This 
absence of fracti onation in the sunspot presumably r elates to the absence of acoustic waves 



in sunspots (e.g. iMug 



strong magnetic field (jParchevsky fe Kosovichevll2007l ). The fact that this is observed in lines 



ach. Hofma nn. fe Staudd 120051) because convection is inhibited by the 



of C I and Fe II suggests that fractionation must set in at lower altitudes than originally 
modeled by Laming (2004; see Figure 1 (left and right panels), where O and C are becoming 
ionized in the region of fractionation, and one would expect neutral O I and C I to be emitted 
from lower, unfractionated layers). 

The existence of fractionation at these low altitudes also offers a possible explanation 
of the observation by Phillips et al. (1994) who found rather small difference in the abun- 
dances of Fe determined from soft X-ray flare plasma, compared with that lower down in the 
atmosphere, determined from the Fe Kf3 fluorescent line, rather than the strong FIP effect 
expected. More recently Murphy & Share (2005) studied 7-ray emission from flares. Pro- 
tons accelerated into the chromosphere by the flare excite 7-ray emission from the ambient 
plasma when its density reaches about 10 14 — 10 15 cm -3 . Element abundances determined 
from the resulting 7-ray spectrum show the presence of a FIP fractionation. The densities 
at which this occurs correspond to the low chromosphere where sound and Alfven speeds 
are approximately equal, and certainly not the Lyman a plateau region where fractionation 
is expected in the Laming (2004) model. A search for FIP fractionation in photospheric 
lines i.e. below the chromosphere (Sheminova & Solanki 1999) reveals very little, if any 
fractionation. Thus all available observational evidence suggests that the low chromosphere 
as another plausible place for FIP fractionation to occur. 

We speculate that the growth of Alfven waves from sound waves near the (3 = 1 layer 
will give an extra ponderomotive force in this region that can account for this. Zaqarashvili 
& Roberts (2006) give a treatment of the parametric conversion of sound waves into Alfven 
waves which requires (3 = 1 when both are traveling in the same direction along the magnetic 
field. This distinguishes it from the phenomenon of mode conversion, which requires nonzero 
wavevector perpendicular to the magnetic field to proceed (e.g. McDougall & Hood 2007), 
and par ametric conversion lower dow n in high f3 plasma where the sound waves must be 



oblique (I Zaqarashvili fe Roberts! 120021 ) . Waves impinging on the chromosphere from below 
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are fast magnetoacoustic waves from the high (3 (gas pressure/magnetic pressure) solar in- 
terior. In the absence of mode conversion, these retain their acoustic character propagating 
as a slow mode wave when (3 « 1 further up in the chromosphere (McDougall & Hood 
2007). At the altitude where (3 ~ 1 (i.e. where the phase speeds of magnetic and acoustic 
waves are similar) these waves can mode convert into other MHD wave modes (Bogdan et 
al. 2003). This would be consistent with the findings of Sheminova & Solanki (1999), who 
find essentially no FIP effect at photospheric altitudes. Acoustic waves will produce no pon- 
deromotive force, and only once mode conversion to the other MHD modes has occurred can 
fractionation proceed. 



Conclusions 



In conclusion then we have refined the model of [Laming (j2004al ) for FIP fractionations 
arising from the ponderomotive force as Alfven waves propagate through the chromosphere. 
We have implemented a non-WKB treatment of the wave transport, which can be further 
modified to include the effects of wave growth and damping, and made a correction to the 
previous formalism to include the Alfven wave transverse velocity in the chromospheric ion 
and neutral partial pressures. The new effects are a saturation of the FIP effect at the correct 
level, and predicted depletions in the coronal abundances of He and Ne, again consistent 
with observations. We find the best match to the observed coronal or solar wind element 
abundances arises for models with an Alfven wave energy fluxes sufficient to heat the corona 
or accelerate the solar wind. The inference that a coronal source of Alfven waves provides a 
FIP effect better matching the observations suggests that coronal abundance anomalies may 
provide novel insights into the coronal heating mechanism(s). 



This work was supported by NASA Contract NNG05HL39I, and by basic research funds 
of the Office of Naval Research. I thank Daniel Savin, Cara Rakowski and an anonymous 
referee for comments on the manuscript. 



A. The Ponderomotive Force 



The ponderomotive force arises from the effects of wave refraction in an inhomogeneous 
plasma. In a nonmagnetic plasma, the refractive index, y/e, is given by e = 1 — uJp/uJ 2 where 
u p is the plasma frequency. Waves are refracted to high refractive index, which means low 
plasma density. The increased wave pressure can then expel even more plasma from the low 
density region, leading to ducting instabilities. In magnetic plasma, e = 1 — u 2 / (u 2 — fl 2 ), 
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where f2 is the ion cyclotron frequency. Thus waves refract to high density regions, and 
plasma is attracted to regions of high wave energy density. A simple expression for the 
ponderomotive force on an ion may be derived as follows. The Lagrangian density for a 
system of thermal plasma of density n with particle mass m and waves is 



L = E \ m * ( V thi,i + VL,) + E 7 ( V «M + V osc,i) • + — — — (Al) 



where is the thermal speed and v osc>i is the oscillatory speed induced by the wave of 
particle i, with mass m«, and charge Wave electric and magnetic fields are given by 
5E and 5B respectively, and 5A is the vector potential. We have omitted the interaction 
term involving the electrostatic potential, since this is constant in a neutral plasma. Putting 
8B 2 /8n = Ei ™v 2 osc j2 + 5E 2 /8n and v OSCji • 5 A = for MHD waves, then 

l = E ?wa+^=^ = E r4,„+E ! V '«' A+ E 2m , ( |-^ 

i i i i i 1 

(A2) 

The 'V Euler-Lagrange equation gives 

d qf d5E 2 

J t (™ th ,, z ) = 2mdnl _ uj2) —, (A3) 

neglecting the spatial variation of B and hence Q iy and evaluating for the component of v thji 
orthogonal to A and B. This is the same as the expression derived by Landau, Lifshitz & 
Pitaevskii (1984), and agrees with earlier work (e.g. Lee & Parks 1983) if 5E 2 = 5E 2 /2, 
where SE P is the peak electric field in the wave, giving a ponderomotive force 



5E 2 



q, 



1 ME * (A4) 



4m t {Q 2 -uj 2 ) dz ' 

When oj « Qi, the ponderomotive acceleration is thus independent of ion mass, which 
is one crucial property relevant to obtaining an almost mass independent fractionation as 
observed. It is also independent of ion change, so long as the ion is charged (and not neutral). 
Litwin & Rosner (1998) give a similar expression derived from the j x B term in the MHD 
momentum equation. 



B. The Non-WKB Transport Equations 

We start from the linearized MHD force and induction equations, 

dSv . . (Vx5B)xB (B • V) SB - (V£B) • B 
P — + V (pu ■ *v) = J = * 4 / , (Bl) 
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and 

= V x (Sv x B) + V x (u x SB) = (B • V) Sv - SBV • u - (u • V) SB, (B2) 

at 

where u and B are the unperturbed velocity and magnetic field, Sv and SB are the pertur- 
bations, and p is the density. Equation (Bl) is rewritten using V (pu • Sv) = pu x V x Sv + 
Sv x V x (pu) + (pu • V) Sv + (5v • V) pu to yield 



dSv 



+ (u ■ V) f v = Vx • V f-gJ) + -f=^ + ( ™> ^ B - <V ' V (B3) 

VV^rpy V^rp 2p 47rp p 



where V A = B/v^ is the Alfven velocity. Writing (VB)-5B = (dB x /dx) SB = - (dB z /dz) SB/2 
since V • B = (assuming dB x /dx = dB y /dy), and similarly for (Vpu) • Sv, and using 

d (pu z /B z ) jdz = gives 



dSv 



, ^xr „/ 5B \ 5B Va SB V a . u 

Here 1/H B = d\nB z /dz, l/H D = dlnp/dz, and below 1/H A = d\nV A /dz. Similar manip- 
ulations give the induction equation in the form 

d ( SB \ , . SB _ „ x r SB u . V A SB u 

Taking equation (B4) plus or minus equation (B5) and rearranging gives the final result, 

9I± + (« ± V A ) ?± = (u ± V A ) + (B6) 



dt v ' dz v ^ \4H D 2H A J ' 



where 7± = Sv ± SB/y/Anp, representing waves propagating in the =F z-directions. 
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Table 1: Coronal FIP Fractionations 



ratio 




relative wave energy flux 






obs. 






1/64 


1/16 


1/4 


1 


4 


16 


64 


a 


b 


c 


He/H 


1.0 


1.1 


1.1 


0.79 


0.55 


0.46 


0.44 


0.68 






C/H 


1.1 


1.3 


1.5 


1.2 


0.84 


0.66 


0.60 


1.36 






N/H 


1.1 


1.3 


1.5 


1.2 


0.86 


0.72 


0.69 


0.72 






O/H 


1.1 


1.4 


1.9 


1.8 


1.3 


1.1 


1.1 


1.00 






Ne/H 


1.1 


1.2 


1.4 


1.1 


0.76 


0.63 


0.60 


0.58 






Na/H 


1.2 


1.9 


4.6 


9.8 


13. 


13. 


13. 




7.8 




Mg/H 


1.2 


1.7 


3.5 


6.1 


6.9 


6.6 


6.3 


2.58 


2.8^ 


3.0+- 


Al/H 


1.2 


1.6 


2.8 


3.6 


3.1 


2.7 


2.5 




3.6+- 


6-8^ 


Si/H 


1.2 


1.7 


3.0 


4.2 


3.8 


3.2 


3.1 


2.49 


5-1-?. 




S/H 


1.1 


1.4 


2.0 


1.9 


1.4 


1.1 


1.0 


1.62 


2.2±0.2 


o 0+1.3 


Ar/H 


1.1 


1.4 


1.8 


1.5 


1.1 


0.91 


0.87 






K/H 


1.3 


2.3 


8.0 


25. 


38. 


41. 


42. 






4 9 +6 - 3 

^ -2.5 


Ca/H 


1.2 


1.6 


2.6 


3.0 


2.3 


1.9 


1.8 






o 1+1.8 


Fe/H 


1.2 


1.6 


2.8 


3.3 


2.6 


2.2 


2.1 


2.28 


4.4±0.5 




Ni/H 


1.2 


1.6 


2.4 


2.5 


1.9 


1.5 


1.5 








Kr/H 


1.1 


1.4 


1.9 


1.7 


1.3 


1.1 


1.0 








Rb/H 


1.3 


2.3 


7.4 


19. 


25. 


25. 


25. 








W/H 


1.3 


2.3 


6.9 


16. 


18. 


17. 


17. 









Note. — FIP fractionations corresponding to the chromospheric model in Figure 11, for relative wave 
energy flux = 1. Other relative wave energy fluxes are given to show the behavior of the FIP fractionation. 
Observational ratios are taken from, (a) Zurbuchen et al. (2002), relative to O, (b) Bryans et al. (2008), 
relative to the mean of O, Ne and Ar, and (c) Giammanco et al. (2008), relative to H. 
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Table 2: Coronal Hole FIP Fractionations 



ratio 


1/04 


relative wave energy 
1/16 1/4 1 4 


nux 

10 


04 


obs. 


rie/li 


l.U 


1.0 


1.0 


1.0 


0.95 


0.92 


0.91 


0.58 


U/H 


l.U 


1.0 


1.1 


1.1 


0.99 


u.yo 


U.yo 


i/ii 
1.41 


i\t /tr 
IN / rl 


l.U 


1.0 


1.1 


1.1 


1.0 


u.y i 


u.yo 


n no 

u.yo 


U/H 


l.U 


1.1 


1.1 


1.2 


1.1 


l.U 


l.U 


i nn 
l.UU 


Mo /FT 


i n 

l.U 


1.0 


1.1 


1.0 


0.98 


n Q/i 


u.yo 


n A7 


Na/H 


1.2 


1.3 


1.4 


1.4 


1.3 


1.2 


1.2 




Mg/H 


1.1 


1.2 


1.3 


1.2 


1.2 


1.1 


1.1 


1.92 


Al/H 


1.1 


1.1 


1.1 


1.1 


1.1 


1.0 


1.0 




Si/H 


1.1 


1.1 


1.1 


1.1 


1.1 


1.0 


1.0 


1.86 


S/H 


1.0 


1.0 


1.1 


1.1 


1.0 


0.98 


0.97 


1.56 


Ar/H 


1.0 


1.1 


1.1 


1.1 


1.0 


1.0 


0.99 




K/H 


1.2 


1.2 


1.3 


1.3 


1.2 


1.2 


1.2 




Ca/H 


1.0 


1.1 


1.1 


1.1 


1.1 


1.0 


1.0 




Fe/H 


1.0 


1.1 


1.1 


1.1 


1.1 


1.0 


1.0 


1.67 


Ni/H 


1.0 


1.1 


1.1 


1.1 


1.1 


1.0 


1.0 




Kr/H 


1.0 


1.1 


1.1 


1.1 


1.1 


1.0 


0.99 




Rb/H 


1.2 


1.2 


1.3 


1.3 


1.3 


1.2 


1.2 




W/H 


1.1 


1.1 


1.2 


1.2 


1.1 


1.1 


1.1 





Note. — FIP fractionations corresponding to the chromospheric model in Figure 15, for relative wave 
energy flux = 1. Other relative wave energy fluxes are given to show the behavior of the FIP fractionation. 
Observational results are taken from Zurbuchen et al (2002), given relative to O. 
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Fig. I. — Force free magnetic field, computed from Athay (1981), from the center of a 
network segment (x = 0) to the center of a supergranule cell (x = 1). We take x = 1 
to represent 1000 km, and y=0 to represent an altitude of 500 km above the photosphere. 
The solid lines represent magnetic lines of force, and dashed lines are logarithmically spaced 
contours of the Alfven speed, assuming the density falls off exponentially with height. FIP 
fractionation in this work occurs towards the top of the chromosphere, where the magnetic 
field is nearly parallel. 
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Fig. 2. — Cartoon illustrating the model. Alfven waves are incident on the coronal loop 
from below on the right hand side. Waves are either transmitted into the loop or reflected 
back down again. Waves in the coronal loop bounce back and forth, with some leakage at 
each footpoint. The magnetic field is taken to be uniform in the coronal section of the loop 
(illustrated), while it varies according to Figure 1 within the chromosphere (not shown on 
this figure). 
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Fig. 3. — Coronal section of loop, length 100,000 km, magnetic field 9.9 G, (half wavelength 
long) showing from top: Elsasser variables in km s" 1 (5B/y/4irp solid lines, 5v dashed lines), 
with black lines for real parts and gray lines for imaginary parts. The loop is approximately 
half a wavelength long. Middle; wave energy fluxes in ergs cmT 2 s _1 , the thin solid line 
shows the difference in energy fluxes divided by the magnetic field strength and should be 
a horizontal line if energy is properly conserved. Bottom, the ponderomotive acceleration 
in cm s~ 2 . Positive acceleration means positive along the z axis, which is upwards pointing 
near z = and downwards near z = 100, 000. 
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Fig. 4. — Same as figure 3 giving the first three panels for the left hand chromosphere 
"A" , where waves leak down from the corona. The extra bottom right panel shows the FIP 
fractionation for Fe, O, and He, relative to H. 
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Fig. 5. — Same as figure 4 for the right hand side chromosphere "B", where Alfven waves 
are launched up from the convection zone. 
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Fig. 6. — Coronal section of loop, length 100,000 km, magnetic field 19.8 G, showing from 
top: Elsasser variables in km s _1 (SB/^Anp solid lines, 8v dashed lines), with black lines 
for real parts and gray lines for imaginary parts. The loop is now a quarter wavelength long 
and reflects most Alfven waves incident from below, and consequently has much smaller 
nonthermal motions than the previous resonant case. Middle; wave energy fluxes in ergs 
cm -2 s _1 . Bottom, the ponderomotive acceleration in cm s~ 2 . 
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Fig. 7. — Same as figure 6 giving the first three panels for the left hand chromosphere 
"A" , where waves leak down from the corona. The extra bottom right panel shows the FIP 
fractionation for Fe, O, and He, relative to H. In the case of a loop off resonance, no waves 
are transmitted through to chromosphere "A" , and no fractionation occurs. 
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Fig. 8. — Same as figure 7 for the right hand side chromosphere "B" , where Alfven waves are 
launched up from the convection zone. Almost complete reflection of Alfven waves occurs 
from the loop footpoint, leading to no fractionation. 
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Fig. 9. — Same as figure 4 (on resonance case) giving the first three panels for the left hand 
chromosphere "A" , where waves leak down from the corona. The extra bottom right panel 
shows the FIP fractionation for Fe, O, and He, relative to H. The wave energy flux has 
been increased by a factor 20, leading to stronger coronal nonthermal motions, and stronger 
fractionation. Helium is now depleted in the corona relative to the chromosphere. 
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Fig. 10. — Same as figure 5 (on resonance case) for the right hand side chromosphere "B", 
where Alfven waves are launched up from the convection zone. The wave energy flux has 
been increased by a factor 20, leading to stronger coronal nonthermal motions, and stronger 
fractionation. 
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Fig. 11. — Same as figure 8 (off resonance case) for the right hand side chromosphere "B", 
where Alfven waves are launched up from the convection zone. The wave energy flux has 
been increased by a factor 20, leading to stronger coronal nonthermal motions, and stronger 
fractionation. Waves are now reflected, and a small inverse FIP effect results. 
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Fig. 12. — Same as figure 9 giving the first three panels for the left hand chromosphere 
"A" , where waves leak down from the corona. The extra bottom right panel shows the FIP 
fractionation for Fe, O, and He, relative to H. Three wave frequencies are now introduced to 
simulate more nearly a realistic chromospheric power spectrum. 
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Fig. 13. — Same as figure 11 for the right hand side chromosphere "B", where three Alfven 
waves are launched up from the convection zone. The FIP fractions are reduced from those 
with a single incident wave, even though one of the waves here is on resonance. The contri- 
butions to the partial pressure of the other waves "dilute" the fractionation, by increasing 
the value of in the denominator of the integral in equation 6, without increasing the 
ponderomotive acceleration. 
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Fig. 14. — Coronal section of open field region up to 500,000 km altitude, showing from 
top: Elsasser variables in km s _1 (SB /y/Anp solid lines, 5v dashed lines), with black lines for 
real parts and gray lines for imaginary parts. Middle; wave energy fluxes in ergs cm -2 s -1 . 
Bottom, the ponderomotive acceleration in cm s~ 2 . 
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Fig. 15. — Same as figure 11 for the right hand side chromosphere "B", where Alfven waves 
are launched up from the convection zone into an open field region. The FIP fractionations 
are evaluated with an incident coronal hole spectrum, as opposed to that for an active region, 
and show the absence of strong fractionation consistent with observations of the fast solar 
wind and coronal holes. 



